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Abstract 

In this work we study the ac dissipation of mesoscopic superconductors at mi- 
crowave frequencies using the time dependent Ginzburg-Landau equations. Our 
numerical simulations show that the ac dissipation is strongly dependent on the 
orientation of the ac magnetic field (h ac ) relative to the dc magnetic field (Hd c ). 
When h ac is parallel to Hdc we observe that each vortex penetration event produces 
a significant supression of the ac losses because the imaginary part of the ac suscep- 
tibility as a function of Hdc increases before the penetration of vortices, and then it 
decreases abruptly after vortices have entered into the sample. In the second case, 
when h ac is perpendicular to Hdc, we observe that the jumps in dissipation occur at 
the same values of Hd c but are much smaller than in the parallel configuration. The 
behavior of the dissipation in the perpendicular configuration is similar to previ- 
ous results obtained in recent microwave experiments using mesoscopic litographed 
squares of Pb (1). 
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In macroscopic samples the ac dissipation is generated by the motion of vor- 
tices that are typically located some micrometers away from the sample sur- 
face, while the dissipation due to surface currents itself can be neglected (0). 
The opposite scenario appears in mesoscopic samples where the contribution 
of the surface currents dominates the dissipation, and vortices play a sec- 
ondary role due to its confinement by the surface currents (j3|). Recently, the 
ac dissipation of an array of mesoscopic square superconductors has been mea- 
sured at microwave frequencies (jl|). The resulting ac dissipation presents a non 
monotonous and oscillating dependence with magnetic field that could be as- 
sociated with the modulation of the vortex entrance events on the surface 
dissipation. 
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In this paper we study the importance of the relative orientation of the ac 
and dc magnetic fields (h ac and H dc ) on the surface dissipation of mesoscopic 
superconductors. Fig. 1(a) represents the parallel configuration, h ac \\ H dc , and 
Fig. 1(b) shows the perpendicular configuration, h ac _L H dc . The doted lines 
show planes, perpendicular to h ac , where the induced ac currents circulate. In 
the present study we neglect the influence of the demagnetization field assum- 
ing that the exact treatment thereof only leads to quantitative corrections, 
while the qualitative behavior remains the same. 
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Fig. 1. The studied magnetic field configurations, (a) h ac \\ H dc and (b) h ac _L H dc . 
Doted lines show planes of circulation of the induced ac currents. 

To study the effects of the anisotropy on the ac response we made numeri- 
cal simulations of the time- dependent Ginzburg-Landau equations in three- 
dimensions. The normalized equations are, see Refs. (y) and (jj), 



~dt 
8A 
~dt 



12 



iV - A) 2 m + (l - T) (|^| 2 - 



1 - T)9?e ^*(-«V - A)V - k 2 (V x B), 



where \I> is the order parameter, A is the vector potential and T is the tempera- 
ture. In the above equations, lengths have been scaled in units of the coherence 
length £(0), times in units of t = Ana n \ 2 L / 'c 2 , A in units of if c2 (0)£(0), \I> in 
units of ^ooiT) = [mc 2 /87re 2 A(T) 2 ] 1 / 2 and temperature in units of T c . 



The boundary conditions for this problem are 



(n^) ± = (V-il) ± ^ = 
V x A(t) = H dc z + [hoc cos(ut)]jl, 

where H dc is oriented in the z direction (see Fig. 1) and fi is the direction of the 
ac magnetic field. When fi = z we obtain the parallel configuration of Fig. 1(a) 
and when (i = x or (i = y we obtain the perpendicular configuration shown 
in Fig. 1(b). In this work we study the behavior of the imaginary part of the 
ac magnetic susceptibility x"(H dc ,uj) = / 27r M(t) sin(a;i). The numerical 
method to solve the three-dimensional set of differential equations is explained 
elsewhere, see Refs. (0) and (0). 
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In Figs. 2(a) and (b) we show numerical results obtained for samples of size 
10A x 10A x 1A, where A is the zero temperature penetration depth. We consid- 
ered here k = 2 and T=0.5. Fig. 2(a) shows the imaginary part of the ac sus- 
ceptibility (x") as a function of H dc for the parallel configuration {h ac \\ H dc ). 
We observe that the dissipation decreases after each vortex penetration event, 
and increases continuously when the number of vortices is constant. Vortices 
enter following the square symmetry of the sample in packages of four vor- 
tices (|5|). However, thermal fluctuations and surface imperfections can allow 
the number of vortices to increase one by one as was shown in Ref. (0). The 
results of Fig. 2(a) are valid for perfect samples at zero temperature. A similar 
qualitative behavior has been observed in Ref. (0) where a two-dimensional 
approach was used. 
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Fig. 2. (a) and (b) are numerical results for the two relative orientations between 
hoc and Hdc- I n ( c ) we show recent experimental results obtained in mesoscopic 
samples (taken from Ref. yj)). 

Fig. 2(b) shows results obtained using the perpendicular configuration (h ac _L 
Hdc)- We observe that the vortex penetration events occur at the same mag- 
netic fields as in Fig. 2(a). This is a consequence of the small value of h ac 
compared with H dc (h ac = 5 x 10~ 4 if c2 <C H dc )- In the linear regime, the mag- 
nitude of the dc magnetic field defines in both cases the condition for vortex 
penetration. Apart from this similarity, the general behavior of x" vs - Hdc i n 
Fig. 2(b) is different to the one in Fig. 2(a). In the perpendicular configura- 
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tion (Fig. 2(b)) we see that vortex entrances produce smaller jumps in x" (Hdc) 
compared with the jumps occurring in the parallel configuration (Fig. 2(a)). 
In Fig. 2(b) we only observe a small decrease in x" an d a change of slope at 
H pi (the i-th penetration field). 

Fig. 2(c) allows a direct comparison between numerical and recent experimen- 
tal results. In Ref. (jl]) the ac dissipation of an array of mesoscopic square 
superconductors was measured at microwave frequencies using the perpendic- 
ular configuration. A simple qualitative comparison shows that the experi- 
mental results are better described by the numerical results of Fig. 2(b) which 
corresponds to the perpendicular configuration. Spectra of Ref. {H) were ac- 
quired at 9.5GHz while the numerical simulations corresponds to a frequency 
of uj = 0.09t; o , where v = l/t Q ~ 100 — 1000GHz in dirty superconductors. 

The reason for the appearance of the anisotropy is clear from Figs. 1(a) and 
1(b). In both configurations, above H Pt i, vortices are oriented in the direction 
of the dc magnetic field. The ac magnetic field induce currents, responsible 
of the ac dissipation, circulate in planes perpendicular to h ac . In the parallel 
configuration (Fig. 1(a)) the ac currents circulates on the faces of the sample 
where the nascent vortices (p) are located. These ac currents generate peri- 
odic compresions and expansions of the vortices that are localized inside the 
sample, but the dissipation is mainly due to the nascent vortices. On the other 
hand, in the perpendicular configuration (Fig. 1(b)), the ac currents circulate 
mainly on the top and bottom faces of the sample where the order parameter 
is depleted due to the vortex cores and not due to nascent vortices. The dissi- 
pation is larger than in the previous case and has a different dependence with 
Hdc- In conclusion, the orientation of the vortices along the Hd c direction and 
the change in orientation of the induced ac currents in both configurations 
generate the anisotropy of the response. 
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